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We demonstrate a new type of quantum mechanical correlation where phase modulators at distant 
locations, acting on the photons of an entangled pair, interfere to determine the apparent depth of 
modulation. When the modulators have the same phase, the modulation depth doubles; when 
oppositely phased, the modulators negate each other. 
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A nonlinear crystal pumped by a monochromatic 
laser may generate time-energy entangled photon pairs 
through the process of spontaneous parametric down con- 
version These pairs may be spatially separated into 
two channels and their frequencies and arrival times at 
distant detectors measured. If Observer A measures the 
frequency of her photon, she knows, by energy conserva- 
tion, the exact frequency of the photon measured by Ob- 
server B. If, instead, she measures the arrival time of her 
photon, she can predict the detection time measured by 
Observer B to within a small window that varies inversely 
with the photon bandwidth Q and may be lengthened 
by dispersion The ability to measure relative time 
and relative frequency, with accuracies not limited by 
time-energy uncertainty, is the hallmark of time-energy 
entanglement [1, H, 0|. 

An important consequence of time-energy entangle- 
ment, first noted by Franson [3,(81] and observed by Bren- 
del et al. (oj , is nonlocal cancellation of dispersion. When 
the photons of an entangled pair are sent through differ- 
ent channels having arbitrary dispersions, the dispersion 
in one channel may be negated by dispersion of the op- 
posite sign in the other channel. This effect results from 
quantum interference and has no classical analog. As al- 
ways, nonlocal effects do not imply transmission of clas- 
sical information. 

In this Letter, we report the first observation of a time- 
frequency analog to nonlocal dispersion cancellation and 
term this effect as nonlocal modulation [lo| . Spontaneous 
down conversion is used to generate entangled photons 
that each have a spectral width of about 280 cm^^. Si- 
nusoidal phase modulators, each operating at a frequency 
of 30 GHz (1 cm~^), are placed in the two beams, labeled 
Channel 1 and Channel 2 (see Fig.[T]). When either beam 
is dispersed by a prism or grating and viewed in the fre- 
quency domain (i.e. as a function of position), this modu- 
lation is completely hidden by the much broader spectral 
width of the photon. To observe the modulation, we cor- 
relate in the frequency (spatial) domain. In the absence 
of modulation, a photon detected at a particular posi- 
tion in Channel 1 coincides with a photon of frequency 
UJ2 — i-^p — in Channel 2, and the correlation in the 
frequency (spatial) domain is therefore a delta function 
(5((jJi +102 — ^p)- When synchronously driven modulators 
are placed in the signal and idler channels, this correla- 
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FIG. 1: (color online). Schematic of nonlocal modulation. 



tion becomes a distribution of discrete sidebands spaced 
by the modulation frequency. What is strange and in- 
teresting is that these distant modulators now act cu- 
mulatively. For example, if the two identical modulators 
have opposite phase, they negate each other and act as if 
neither modulator were present. Conversely, if operated 
with the same phase, they produce the same correlation 
as does a single modulator with twice the modulation 
depth acting on only one of the photons. 

To avoid confusion, we mention a different type of 
quantum interference discovered by Steinberg and col- 
leagues ll|, |l2| that bears on Hong-Ou-Mandel interfer- 
ometry. Here, because it is not possible to determine 
which photon passed through a dispersive medium, there 
is an interference of Feynman paths, and even-order dis- 
persive terms are not seen by the interferometer. Recog- 
nizing the importance of time-frequency duality, Tsang 
and Psaltis have suggested the equivalent of the Stein- 
berg interference in the time domain [3, . Since the 
photons meet on a single beam splitter in these examples, 
the interference is a local effect, and classical analogies 
have been demonstrated 
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A schematic of the experiment is shown in Fig. [21 
We pump a 20 mm long, periodically-poled, magnesium 
oxide-doped stoichiometric lithium tantalate crystal (PP- 
SLT, HC Photonics Corp.) with 0.8 W from a 532 nm 
cw laser (Coherent Verdi VIO). The nonlinear crystal is 
phase matched to produce 32 nm bandwidth, degenerate 
photon pairs at 1064 nm. All fields are polarized along 
the extraordinary axis of the crystal. The generated pho- 
tons are filtered from the strong 532 nm pumping beam 
using a four-prism setup and are then coupled into a 
polarization-maintaining fused-fiber beam splitter which 
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FIG. 2: (color online). Experimental setup. and $2(i) 

are 30 GHz sinusoidal modulators synchronously driven with 
variable relative phase (see text for details). 



diverts the photons into Channels 1 and 2 with equal 
probabihty. The photons pass through identical sinu- 
soidal phase modulators (EOSPACE) driven at 30 GHz 
with modulation depths of about 1.5 radians. The rela- 
tive phase between the modulators is controlled using a 
calibrated phase trimmer. Following the modulators are 
identical monochromators, each having a linear disper- 
sion of 210 GHz/mm and a Gaussian instrument response 
function with a FWHM bandwidth of 8.5 GHz. To ob- 
tain frequency-domain correlations, we fix the output slit 
in Channel 1 at Xi and scan the position X2 in Channel 
2. The photons transmitted through the monochromator 
slits are coupled into multimode fibers and detected with 
single photon counting modules (SPCMs, id Quantique 
id400 and PerkinElmer SPCM-AQR-16-FC). 

The primary experimental results of this work are 
shown in Figs. [3] and [H For each case, we set the 
monochromator slit in Channel 1 at an arbitrary posi- 
tion xi near the center of the generated 32 nm spectrum 
and leave the position of this slit fixed thereafter. The 
slit in Channel 2 is scanned over positions X2, and the 
coincidence rate of the two detectors (with gate width 
T — 1.25 ns) is recorded as a function of this position. 
For each position, the rate is averaged for 20 s. 

With the pump frequency defined as LUp, and the po- 
sition X2 proportional to the frequency uj2, we express 
the coincidence rate as a function of relative frequency 
A = 0^2 — {ujp — LUi). The scale of the frequency axis is 
calibrated by measuring the sideband spacing of a single- 
mode 1064 nm laser modulated at 30 GHz, with the zero 
position chosen (at the start of the experiment) as the 
location of the correlation peak for unmodulated photon 
pairs. 

Figure [3] shows the experimental results without mod- 
ulation and with modulation in a single channel. In 
Fig. EJa), both modulators are turned off by discon- 
necting their 30 GHz drive signals. As expected by 
energy conservation, a single correlation peak is ob- 
served. In Figure jS^b), Channel 1 is phase modulated as 



exp[i5cos(wmf)] with a modulation depth of 5 = 1.5, and 
Channel 2 is not modulated. The frequency correlation 
is now distributed over a set of sidebands, having Bessel 
function amplitudes Jni^), whose total area is equal to 
that of Fig. Eta). 
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FIG. 3: Frequency correlation measurements (a) with both 
modulators turned off and (b) with the modulator in Channel 
1 running at a modulation depth of 1.5. Dots are data; curves 
are theoretical fits (see text). All data is approximately shot- 
noise- limited. 



T 


200 










0) 


150 








0) 

(J 


100 


c 




0) 




T3 


50 


'o 


C 




Coi 






(a) 



120 -60 60 120 
A (GHz) 



200 
150 
100 
50 




1 







-120 -60 60 120 
A (GHz) 



FIG. 4: As in Fig. |3] but both modulators running (a) with 
the same phase and (b) with opposite phase. 

In Fig. m^a) , both modulators are turned on at a mod- 
ulation depth of 5 — 1.5, and the cable length is adjusted 
so that they have the same phase. They now act cu- 
mulatively (constructively interfere) to produce a set of 
sidebands having a Bessel function distribution J,^(2(5). 
The frequency-domain correlation function of two distant 
modulators is therefore the same as that which would be 
obtained by correlating an unmodulated photon with a 
photon modulated at twice the modulation depth. 

In Fig. m^b) , the modulators are run at the same depth 
as in the previous paragraph, but now the relative cable 
length is adjusted so that the modulators are run in phase 
opposition. The modulators now destructively interfere, 
and no sidebands are visible. The solid curves in Figs.|3| 
and m are theoretical fits to the data. 

The theory of nonlocal modulation as measured in the 
frequency domain has been developed by Harris [13] for 
the case of frequency correlation using ideal detectors 
with perfect frequency (spatial) resolution. In the fol- 
lowing paragraphs, we develop the theory to allow for 
finite-resolution monochromators and detectors having 
both arbitrary transmission functions and specified tem- 
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poral gate widths. Working in the Heisenberg picture, a 
nonUnear crystal of length L is pumped by a monochro- 
matic laser at frequency ujp. A positive- frequency field 
operator a(a;, z), representing entangled photons, evolves 
inside the crystal and may be written in terms of an en- 
velope 6(ci;, z) which varies slowly along the propagation 
direction: a{u!,z) = 6(a;, z) exp[ifc(ci;)z]. The propagation 
equations describing entangled photon generation are 



db{uj, z) 

dz 
dh\uj,z) 
dz 



— iK(a;)6^(u;p — w, z) exp [iAfc(a;)z], 

= — iK*(cL;)6(wp — cj, z) exp [— iA/c(ci;)z]. (1) 



where k{uj) and Afc(w) are the coupling factor and wave- 
vector mismatch, respectively. The solution for the out- 
put field at z = i, expressed in terms of the vacuum field 
o.\a.c{^) at the input of the crystal, is 



flout (t^) = A(w)avac(w) + B{uj)al^^{ujp ~ u), (2) 

where, to preserve the commutation relations, the func- 
tions A{uj) and B{uj) satisfy - \B{ijj)\^ = 1 and 
A{uj)B{Ljp — Lo) = B{ijj)A{ujp — uj). 

The time-domain output field operator is related to 
its frequency-domain counterpart [Eq. ^] by the inverse 
Fourier transform, aout(i) — /^o^ lout(^) exp(— ia;t)dw, 
and is normalized so that the total rate of generated 
photons exiting the crystal is i?out = (fllut(^)flout(i))- 
The generated photons are separated into two chan- 
nels, denoted as Channel 1 and Channel 2, using a 
50/50 beam splitter. The field operators at the out- 
puts of the beam splitter are ai{t) = 02 (t) = -^aout(i)- 
The photons are modulated by periodic phase modu- 
lators whose time-domain, Fourier-series transfer func- 
tions are mi{t) — '^f,qkexp(—ikLUmt) in Channel 1 and 
"^2(0 — '^irieyi'p{~ilLUjnt) in Channel 2, with Fourier 
transforms mi{(jj) = J2k1k^i^ ~ kuim) and m2{w) — 
J2i riS{uj — lujm), respectively. With the * symbol denot- 
ing convolution, the frequency-domain modulated fields 
are ai{Lu) = ai{uj) * mi(a;) and 02(0;) = 02 (w) * 1712(10). 
Substituting ai{uj), a2(w), mi(a;), and m2(w) into the 
expressions for ai(aj) and 0,2 (lo) yields 



ai{uj) 



02 (w) 



^ 00 

^ qk[A{uj - fcwm)avac(t^ - kuj„i) 

k— — oo 

+ B{uj ~ ku}yn)a\,^c{^p - w -f kujjn)] , 
^ 00 

— = ^ n [A{uJ - luJm)avac{uJ ~ luJm) 
^ l=~oo 

+ B{UJ - luJm)al^^{uJp - UJ + iLOyn)] ■ 



The modulated photons are frequency correlated by 
passing each through identical monochromators whose 
output slits may be translated to select frequencies uoi — 
(3xi in Channel 1 and uj2 = (3x2 in Channel 2, where 
the constant f3 is the linear dispersion of the grating 
systems. The monochromators (spectral filters) have 
field transmission functions Hi(lu — f3xi) and H2{uj — 
(3x2)- The filtered field operators in Channels 1 and 
2 are aii{uj,xi) — ai{uj)Hi{uj — l3xi) and a2i{uj,X2) = 
d2{uj)H2{uj — (3x2), respectively. The count rates at the 
outputs of the monochromators are given by i?i(a;i) = 
{a\f{t,xi)aif{t,xi)) and -R2(a;2) = {a\i{t,X2)a2i{t,X2)). 
These rates are 
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X / \B{uj — kujm)\^ \Hi{uj — f3xi)\^d{jj, 
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X / \B{UJ - luJrn)\^ \H2{UJ - pX2)\^d0J. (4) 



Assuming a gate width T, the coincidence rate 
for the two detectors is related to the second- 
order Glauber correlation function G^^\ti,xi,t2,X2) = 



{alf{t2,X2)alf{ti,Xi)aif{ti,Xi)d2i{t2,X2))- With the as- 
sumption that the resolution of the monochromators is 
high, or equivalently that the filter widths are small (as 
compared to the modulation frequency Wm), it can be 
shown that the correlation function depends only on the 
difference of the arrival times t — t2 — ti, and the coin- 
cidence rate is 



rT/2 

Rc{xi,X2)= G^'^\T,Xi,X2)dT. (5) 

J~T/2 

Equation ([5|) may be expanded using Wick's theorem 
and shown to be given by 



Rc{xi,X2) = Rl{xi)R2{x2)T 

OC 

y^^qkrn~kFk{T, xi,X2) 



k—-—oo 



dT, (6) 



where A — f3{xi + X2) — tOp, n — [A/uJm + 5J , and 



Fk(T, Xi,X2) 



1 



A(lU — kuJm)B{ujp — + kuJm) 



(3) 



47r Jo 

X Hi{uj ~ l3xi)H2{u}p - UJ ~ (3x2 + nuj,n) 
X exp{iujT)duj. (7) 
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The first term in Eq. ([6]) is the result of accidental 
coincidences between unpaired photons in a gate width T. 
The second term is the coincidence rate between paired 
photons and captures the modulation effects described 
in this Letter. To obtain Eqs. (H])-®, we have assumed 
that the transmission widths of the monochromators are 
small as compared to the modulation frequency and large 
as compared to the inverse of the temporal gate width 
T. In our experiment these assumptions are satisfied by 
factors of 3.5 and 11, respectively. 

If we further assume that A(uj) and B{ujp — uj) are 
constant in the vicinity (±150 GHz in Figs. [3] and [4]) of 
Lu = j3xi and are equal to Aq and i?o, respectively, then 
Eq. (ID) becomes 



i?c(A) = R1R2T + CnH^^\nuJ„r - A), 

where H^^'){uj) = \Hi{lu)\^ * l^2(t^)P, and 
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The solid curves in Figs. [3] and [4] are theoretical fits to 
the data using Eq. ([S]) shifted horizontally so as to match 
center. The Fourier series coefficients for sinusoidal phase 
modulators are Bessel functions with qk = Jfc(— (5i) and 
fi = Ji{—52), where 5i and 82 are the modulation depths 
in Channels 1 and 2, respectively {\5i\ = \52\ = 1.5 
in our experiment). We model the monochromator re- 
sponse functions in Channels 1 and 2 as Gaussians with 
FWHM bandwidths F: Hi{uj) = ai exp [-2\n{2)uj^ /T^] 
and H2{ijj) = a2 exp 

[-21n(2)wVr^]- (The monochro- 
mator in Channel 1 is the mirror image of the one in 
Channel 2 which has a measured FWHM bandwidth of 
8.5 GHz.) The transfer functions include fitting parame- 
ters ai and a2 used in Figs. [3] and |4] to account for trans- 
mission losses and the difference in detection efficiencies 
of the photon counters. 

To obtain the constants ^0 and Eq, for each case in 
Figs. [3]and[4l we measure the average value of R2 and 
use Eq. (j8bp to calculate |-Bol- We obtain j^ol from the 
commutator-preserving condition |^oP ~ l-^oP = 1- For 
all curves, the fitting parameters are taken as a\ — 1.20 x 
10~^ and ^2 = 5.59 x 10~*. These values are in good 
agreement with loss measurements and estimates of the 



photon counter detection efficiency, where we note that 
the id400 detector in Channel 1 has a detection efficiency 
an order of magnitude larger than the SPCM-AQR-16- 
FC detector in Channel 2. 

In summary, this work reports the first observation 
of a quantum effect termed as nonlocal modulation. 
We have experimentally shown how distant modulators, 
when correlated in the frequency domain, may inter- 
fere constructively or destructively. Though this work 
has dealt with the effects of synchronously-driven sinu- 
soidal modulators, a more general statement for nonlocal 
modulation is that phase modulation in Channel 1 of 
the form exp[i$i(i)] acts cumulatively with modulation 
exp[j$2(0] in Channel 2 so as to produce a frequency do- 
main correlation proportional to the square of the Fourier 
transform of exp{z[$i(i) H- <I>2(i)]}. For this relation to 
hold, it is required that <&i(i) and $2(i) both vary slowly 
as compared to the temporal width of the biphoton wave 
function. 
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